The influence of strong magnetic fields on proto-quark stars 
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We use the MIT bag model to analyze different stages of magnetized quark star evolution in- 
corporating baryon number conservation and an anisotropic energy momentum tensor. The first 
stages of the evolution are simulated through the inclusion of trapped neutrinos and fixed entropy 
per particle, while in the last stage the star is taken to be deleptonized and cold. We find that, 
although strong magnetic fields allow for more massive quark stars, the evolution of isolated stars 
needs to be constrained by fixed baryon number, which lowers the star masses. Moreover, magnetic 
field effects, measured by the difference between the parallel and perpendicular pressures, are more 
pronounced in the beginning of the star evolution when there is a larger number of charged leptons 
and up quarks. We also show that having a spatially varying magnetic field allows for larger mag- 
netic fields to be supported since the model employed generates large magnetic fields only at high 
densities, where the longitudinal matter pressure is large enough to partially compensate for the 
negative magnetic field longitudinal pressure. 

PACS numbers: 26.60.kp, 26.60,+c, 95.30.Tg, 95.85.Sz 
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I. INTRODUCTION 

Neutron stars are compact objects with masses of the 
order of 1-2 Mq, radii around 10 km, and a temperature 
of approximately 10 11 K at birth, cooling rapidly by emit- 
ting neutrinos. In conventional neutron star models, they 
are comprised of hadronic matter and leptons. Never- 
theless, according to the Bodmer-Witten conjecture, the 
true ground state of matter could be deconfined quarks in 
environments where densities are very large [l|-|3(. This 
hypothesis lead to the possibility that neutron stars could 
also be strange stars or quark stars, comprised of quarks 
and leptons [4j. The true nature of these compact stars 
remains a source of speculation and further constraints 
are necessary so that models can be ruled out. 

Magnetars have extremely high magnetic fields of up 
to B = 10 15 G on the surface @-@|. Nevertheless, all 
pulsars have relatively strong magnetic fields [{J [To[ , and 
for this reason, any complete analysis of pulsar features 
should include magnetic field effects. Unfortunately, all 
realistic calculations of magnetic field limits in the center 
of stars are model dependent to some extent, even when 
the virial theorem is considered. Some results along this 
line can be found in Rcfs. and have limits rang- 

ing between B = 10 17 — 10 20 G. At such high magnetic 
fields, the thermodynamical as well as the hydrodynami- 
cal properties of matter become anisotropic [19l425l | and, 
therefore, need to be carefully handled. 

Many pulsars exist in binary systems, however, it is 
possible to have an isolated pulsar, in which case one 
needs to take care to properly enforce baryon number 
conservation. This happens because, differently from pul- 
sars in binary systems which can accrete matter from 
a companion, isolated stars cannot and, therefore, they 
conserve baryon number. In this case, as was previously 
shown in Refs. (26l - [28| . there exist windows of stability 



which only allow for stars with particular values of baryon 
number to exist at any given time. The different snap- 
shots of the evolution are simulated though different en- 
tropies per particle and trapped neutrinos. As was shown 
in Ref. [29j for a hadronic model, such an analysis shows 
that the maximum mass that can possibly be attained 
by hadronic stars decreases with time. 

There have been some real-time calculations of the 
evolution of proto-neutron stars performed in Refs. [30l — 
|33| utilizing Boltzmann transport equations with realistic 
neutrino cross sections. Although a very interesting pre- 
liminary calculation of the neutrino absorption cross sec- 
tion influenced by magnetic field effects was performed in 
Ref. [lil for neutron stars, to the best of our knowledge, 
such early time simulations have not been performed in- 
cluding magnetic field effects in quark matter. 

There have been previous studies of proto-quark and 
proto-hybrid stars [3514471 . Others have considered proto- 
neutron stars subject to a magnetic field H?£ |49| and 
zero temperature quar k stars in a magnetic field using 
the bag model [50l l5l|. In this paper we make the first 
study of the effects of magnetic fields on finite tempera- 
ture proto-quark stars incorporating the effect of baryon 
number conservation. This is important when consider- 
ing isolated proto-quark stars. In addition, we include 
the pressure anisotropies generated by all relevant de- 
grees of freedom (quarks, leptons, and neutrinos) and 
the magnetic field. 

Using the MIT bag model [52[ , we determine the effect 
of strong magnetic fields on three snapshots of quark star 
evolution. Our task includes a comprehensive investiga- 
tion of the pressure anisotropy at zero temperature and 
at fixed entropy per baryon, including finite temperature 
effects. In most of the previous works on neutron stars 
subject to strong magnetic fields, the anisotropic pres- 
sures were either disregarded j53l457| or only considered 
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at zero temperature [20, [2l|, [23| . Here we ignore the ef- 
fects of the quark anomalous magnetic moment (AMM) 
and focus on the effect of Landau quantization only. 

The structure of the paper is as follows. In Sec. [Hi we 
review the formal setup of the problem and present ex- 
pressions for the bulk properties of the system. In Scc. lIIII 
we present our results for the transverse and longitudinal 
pressures and analyze the effect on the mass-radius rela- 
tionship and star properties at three different moments 
in its evolution. In Sec. IIVI we present our conclusions 
and an outlook for the future. 



II. FORMALISM 

The Lagrangian density for the bag model including 
the magnetic field is given by 

C = [*, (i 7 % - e q ^A„ - m q ) * 9 - B~\ Q v 

- ^-F^F^ + *i {i-fd M - erfA M - m,) (1) 

107T 

where ty q stands for the 3 light quark fields, e = 0.08 is 
the electron charge (in Gaussian natural units) , A 1 * is the 
electromagnetic field and 0y is a step function defined 
as 1 inside the bag and outside. The bag constant 
value is chosen to be B = (154 MeV) 4 in accordance with 
an analysis of stability made for magnetized proto-quark 
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field tensor (multiplied by 1/V47T due to our choice of 
units). Finally, stands for the electron, muon, and 
electron and muon neutrino fields. 

The magnetic field is taken to be parallel to the z direc- 
tion. For details of the calculation of the bulk prope rties 
of a Fermi gas we refer to Rcfs. [H M, HI El, MM, M~ 
l62j . The pressure and energy density in the bag model 
are the pressure and energy density of a non-interacting 
Fermi gas with the bag constant subtracted and added, 
respectively. For the charged particles, the matter con- 
tribution to the pressure and energy density for finite and 
zero temperature in the presence of a constant magnetic 
field, B, in the z direction are 
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where 7^ is the degeneracy of each particle i (6 for quarks, 
2 for electrons and muons and 1 for neutrinos) taking 
into account spin and/or number of colors, Q ei is the 
electric charge of each particle (in multiples of the elec- 
tron charge), B the magnetic field, k Fi = [i\ — mf 
the Fermi momentum of each particle modified by the 
magnetic field, ifii = \/mf + 2\Q ei \eBv the mass of 
each particle modified by the magnetic field with 
being the bare mass of each particle (with m u ^ — 5 
MeV and m s = 150 MeV). The distribution functions 
for the particles (/+) and anti-particles (/_) are f± t = 
l/{ e W^+™lT^)/T + 1 j j whcre ji is thc temperature. 

The chemical potential of each particle Mi is calcu- 
lated from thc independent chemical potentials (baryon, 
charge and lepton chemical potentials), each of which is 
related to a conserved quantity (baryon number, electric 
charge and lepton fraction) 



M; = QbjMb + Qe t {^q + Mi) + QliHU 



(0) 



where Qsi is the baryon number (1/3 for quarks and 
for leptons) and Qi t thc lepton number (0 for quarks 
and 1 for leptons) of each particle. Note that the lep- 
ton chemical potential is non-zero only when the lepton 
fraction is fixed. If strangeness were conserved, like in 
the case of the initial stages of a heavy ion collision, a 
strange chemical potential would have to be introduced. 

Due to energy level quantization in directions perpen- 
dicular to thc magnetic field, the thermodynamic quan- 
tities for charged particles are summed over the Landau 
levels v, which for zero temperature have a maximum 



value equal to v„ 



L(m^ -mj)/(2\Q ei \eB)\. Note 



that the degeneracy of the zeroth Landau level is always 
half of the usual value for spin 1/2 particles. For non- 
charged particles, such as neutrinos, the pressure and 
energy density expressions take the usual form 52 1. 

For the star calculations, we require that the matter 
is charge neutral and in chemical equilibrium. The first 
condition can be written as 



QeiPi = 0, 
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where the density of each charged particle for finite and 
zero temperature cases is 

P* = E ^2 \Q*i\ eB f (f+i ~ f-i)&> ( 8 ) 
Pi {T = Q) = Y,^\Qei\eBk Fh (9) 



and the total baryon density is simply ^2iQi3iPi- The 
second condition, corresponding to chemical equilibrium, 



can be written by replacing the baryon number, electric 
charge, and lcpton number for each particle in Eq. (|6|): 
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(10) 

It is easy to see from Eq. (|10p that = /i s = fj, u -\-fj, e —fii, 
He = (ip, and \x Vc = \i v ^. 

Due to the magnetic field, the perpendicular compo- 
nent of the matter contribution to the pressure is also 
modified. In fact, one can show that the difference of the 
parallel and perpendicular components of the pressure is 
proportional to the magnetization. The magnetization 
can be calculated by taking the negative of the deriva- 
tive of the grand potential with respect to the magnetic 
field, resulting for finite and zero temperature in 
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so the perpendicular component of the matter contribu- 
tion to the pressure becomes 
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Due to the presence of the magnetic field, there are, 
in addition, contributions from the electromagnetic field 
tensor (see Eq. (Q])) to the pressures and energy density 
[63l HH . Due to the breaking of rotational symmetry by 
the magnetic field, these contributions are again different 
in the parallel and perpendicular directions. The result- 
ing pressures and energy density are 
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where the Air comes from the choice of Gaussian natural 
units. 

The temperature is not expected to be constant in the 
interior of compact stars Q. This can be easily taken 
into account by computing quantities at fixed entropy 
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FIG. 1. (Color online) Binding energy versus baryon density 
for charge neutral matter at zero temperature shown with and 
without (constant) magnetic field effects. For the magnetized 
case B = 4.30 x 10 ls G. 



per particle. Numerical simulations [301 - 13 21 .1 651] show that 
this ratio can reach values 1 or 2, depending on the stage 
of evolution. In this case, the temperature increases nat- 
urally from a smaller value on the surface to a higher 
value in the center of the star. The entropy per particle 
(baryon) can be calculated through the thermodynamical 
expression 

— = — = e + p - ^bPb ( 7) 

A p B Tp B [ > 

Immediately after the supernova explosion, the elec- 
tron and muon neutrinos are still trapped, which can be 
modeled using a fixed lepton fraction defined as 

Ei QhPi 



Y, = 



PB 



which according to simulations |3C 
assume Yj = 0.4 throughout. 



(18) 

is Yi ~ 0.4. Here we 



III. RESULTS AND DISCUSSIONS 

We begin our analysis with the binding energy for mat- 
ter without leptons. The results at zero temperature are 
shown in Fig. [T] In this work, as previously stated, we 
use a bag constant of B = (154 MeV) 4 , which repro- 
duces a minimum in the binding energy {B/A = 925.9 
MeV) slightly below the one for the iron. We point out 
that the inclusion of the magnetic field makes matter 
even more stable, giving a binding energy that is even 
lower (B/A = 914.30 MeV). Such a behavior was already 
pointed out in Refs. jEUIEzj]- 

We now proceed with our analysis of the equation of 
state of the matter only (the pure magnetic field contri- 
butions of ±B 2 /8ir are not included at first). The paral- 
lel and perpendicular pressures of matter are plotted in 
Fig. [5] as a function of energy density of matter for three 
different cases: 




FIG. 2. (Color online) Different matter contributions to the 
pressure versus matter contribution to the energy density at 
different snapshots of the star evolution shown with and with- 
out (constant) magnetic field effects. For the magnetized case 
B = 4.30 x 10 18 G. 



i) s/p B = l,Y t = 0.4, 

ii) s/p B = 2, p Ul = 0, 

iii) s/pB = 0, p Ul = 0. 

These correspond to three snapshots of the time evolu- 
tion of a quark star in its first minutes of life. Such 
an analysis is important, as it has been shown in cal- 
culations/simulations that quark matter can already be 
formed during supernova explosions [(38|,[69[. At first, the 
star is relatively warm (represented by fixed entropy per 
particle) and has a large number of trapped neutrinos 
(represented by fixed lepton fraction). As the trapped 
neutrinos diffuse, they heat up the star [27j. Finally, the 
star can be considered cold. 

Notice in Fig. [2] that (constant) magnetic field effects, 
seen by the difference between the parallel and perpen- 
dicular pressures, are more pronounced in the beginning 
of the star evolution (i), when there is a larger number of 
negatively charged leptons and up quarks. Other effects, 
such as the Haas- van Alphen oscillations (70l - [73| can only 
be seen at the last stage of evolution (iii) , since such os- 
cillations are smoothed out by thermal effects. These 
oscillations can be better seen in Fig. [3] where the par- 
allel and perpendicular pressures of matter are plotted 
as a function of magnetic field for different chemical po- 
tentials representing the center of the star (pb = 1500 
MeV) and an outer region (pb = 1200 MeV). Our results 
for T = can be compared with the ones presented in 
Fig. 1 of Ref. [23j]. The models are not the same, but the 
same physics is obtained for pb = 1500 MeV (which cor- 
responds to p q = 500 MeV in Ref. [H]). At these high 
densities, the matter contribution to the perpendicular 
pressure only deviates significantly from the parallel one 
when the magnetic fields are as high as 3 x 10 19 G. 

The magnetic field configuration considered in Figs. [1] 
and[^]was constant in the z direction with B = 4.30 x 10 18 
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FIG. 3. (Color online) Parallel (red/orange lines) and perpen- 
dicular (respective green/dark green lines) pressures of matter 
versus magnetic field representing different snapshots of the 
star evolution. 



G. One expects a non-constant magnetic field whose mag- 
nitude decreases as one goes out from the center towards 
the edge of the star. To allow for this, we consider a 
magnetic field with a variable strength B*(ps) that in- 
creases with chemical potential (density). It ranges from 
a lower value set to B sur f = 10 15 G at small baryon 
chemical potentials (densities) to a higher value set to 
B c = 4.30 x 10 18 G at high baryon chemical potentials 
(densities) following the prescription from Refs. 



B*(pb) = B sur f + B c 



l-e° 



(19) 



with a = 2.5, b = -4.08 x 10 -4 and p B given in MeV. 
The variable magnetic field prescription employed only 
significantly modifies the behavior of the system at high 
densities and therefore does not change the minimum of 
the binding energy. As explained in [171 ]. stronger mag- 
netic fields are necessary at larger densities in order to 
conserve the magnetic flux. For an alternative explana- 
tion of how the magnetic field strength increases toward 
the center of the star see [77j ■ 

The effects of considering a variable magnetic field can 
be seen in Fig [4] where the EOS is shown. They only 
become visible at high densities. It is important to no- 
tice that the magnetic field in the center of the star never 
reaches B c , but only ~ 80 — 90% B c . This happens be- 
cause the stars become gravitationally unstable for higher 
central chemical potentials. 

We can ask ourselves if it would make sense to increase 
the magnetic field further than we have already done. 
The answer is probably not. Beyond B™ ax = 4.39 x 10 18 
G the parallel pressure (now including the field contri- 
bution) begins to become negative in the center of the 
star, even using the variable magnetic field prescription. 
According to Ref. [78[ this defines an upper limit for the 
magnetic field, beyond which the system becomes unsta- 
ble. Using a constant magnetic field, the upper limit is 
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FIG. 4. (Color online) Different matter contributions to 
the pressure versus matter contribution to the energy den- 
sity at different snapshots of the star evolution shown for 
variable (chemical potential dependent) magnetic field with 
B c = 4.30 x 10 18 G. 



lower, as indicated in Fig. \5\ For a chemical potential 
of fis = 1500 McV the parallel pressure becomes nega- 
tive for B = 3.17 x 10 18 G but for [i B = 1200 McV this 
happens for a lower magnetic field of B = 1.77 x 10 18 
G. Those numbers were obtained for stage (iii) of the 
star evolution. In the earlier stages, the parallel pres- 
sure becomes negative for slightly higher magnetic fields, 
B = 3.37 x 10 18 G for (i) and (ii) at /ig = 1500 MeV. 
As this demonstrates, having a variable magnetic field 
prescription improves the stability since it produces high 
magnetic fields only at high densities, where the longitu- 
dinal matter pressure is large enough to partially com- 
pensate for the negative pure magnetic longitudinal pres- 
sure. 

Having obtained the EoS for the system including mag- 
netic field effects we can use this as input to the TOV 
equations [zi,^ as long asB< B^**. We will use the 
transverse pressure to solve TOV equation. This should 
be a good approximation when P± ~ Py ; otherwise, one 
should also include the breaking of spherical symmetry in 
the solution of the gravitational metric. Wc will consider 
two different central magnetic fields of B c = 2 x 10 18 G 
and B c = 4.3 x 10 18 G, the latter of which is close to 
B™ ra . The results can be seen in the mass-radius dia- 
gram shown in Fig. [6l Without magnetic field effects, 
the maximum mass decreases with time from (i) to (ii) 
as the star deleptonizes and becomes more isospin sym- 
metric (due to charge neutrality). Actually, the star has 
more up than down quarks at (i). From (ii) to (iii) the 
star's maximum mass decreases again as it becomes cold 
and loses thermal energy. On the other hand, when the 
(variable) magnetic field effects are included, the maxi- 
mum mass increases with time, from (i) to (ii) and from 
(ii) to (iii) . The difference in the possible star masses be- 
tween different snapshots is also larger and all stars are 
more massive. The bulk of these effects come from the 
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FIG. 5. (Color online) Parallel (red/orange lines) and per- 
pendicular (respective green/dark green lines) pressures ver- 
sus magnetic field representing different snapshots of the star 
evolution. Pure magnetic field contribution B 2 /8tv included 



inclusion of the magnetic field in the pressure and energy 
density. 

The analysis above works only for stars that are in bi- 
nary systems, and therefore can change baryon number 
with time. In order to analyze isolated stars, we need to 
constrain the star baryon number by the one at its first 
moment of evolution (i). As can be seen in Table [T] the 
maximum masses of stars in all cases now decrease with 
time. This result agrees with similar hadronic calcula- 
tions [29j]. The corresponding radii, on the other hand 
increase with time. It is interesting to notice that within 
this analysis, the magnetic field in the center of stars from 
stages (ii) and (iii) only reaches about ~ 40 — 50% B c . 
The data from Table U with fixed baryon number is also 
represented by symbols in Fig. [6] 

It is important to note that we assume that the mag- 
netic field is constant in time in the analysis of different 
snapshots of the quark star evolution. Such an assump- 
tion is reasonable considering that the timescale for the 
ma gne tic field decays is of the order of 10 5 — 10 7 years 
[8l| - l84l ] . This should be contrasted with the timescale for 
the cooling of proto-neutron stars which is on the order 
of minutes. 



IV. CONCLUSIONS AND OUTLOOK 

In this paper we have examined the effect of strong 
magnetic fields on proto-quark stars using the MIT bag 
model. Three situations were investigated: (i) fixed en- 
tropy per baryon s/ps = 1 with trapped neutrinos, (ii) 
s /pb = 2 without neutrinos and (iii) a zero tempera- 
ture system with no neutrinos. The pressure anisotropy 
inherent to Fermi gases in external magnetic fields was 
obtained both for fixed magnetic fields and for a scenario 
where the magnetic field varies with the baryon chemical 
potential. The conditions for chemical equilibrium be- 
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FIG. 6. (Color online) Mass-radius diagram for star families 
representing different snapshots of the star evolution shown 
without and with (variable) magnetic field. The dots repre- 
sent the time evolution of the most massive stars considering 
fixed baryon number. The lower part of the figure is not 
shown as it exhibits the usual behavior. 



B = B =£0 

stages i ii iii i ii iii 

M max (M Q ) 1.66 1.65 1.65 1.84 1.85 1.88 

M max (M ) fixed A 1.66 1.58 1.55 1.84 1.73 1.70 
R (km) fixed A 9.19 9.53 9.54 9.14 9.62 9.73 

TABLE I. Maximum masses with and without fixing baryon 
number and corresponding radii for different snapshots of the 
star's evolution with and without the inclusion of a spatially 
varying magnetic field. For the case with magnetic field we 
used B c = 2.0 x 10 18 and B c = 4.3 x 10 18 G. 



tween the quarks, leptons, and neutrinos were enforced. 
In the case of isolated proto-quark stars we incorporated 
the effect of baryon number conservation during its evo- 
lution. 

We demonstrated that strong magnetic fields allow for 
more massive quark stars. However, the evolution of iso- 
lated stars needs to be constrained by fixed baryon num- 
ber, which lowers the star masses. For our final mass 
vs radius results we studied two different central mag- 
netic fields: B c = 2.0 x 10 18 G and B c = 4.3 x 10 18 G. 
In both cases the maximum magnetic field generated in 
the star's interior at stage (i) of the evolution is on the 



order of -B max < 0.8£?™ ax . For the lower value of the 
magnetic field considered, this places us well below B™ ax 
where the longitudinal pressure of the system becomes 
negative. In addition, the level of pressure anisotropy in 
this case is relatively small with Pi\/P± ~ 0.85, giving us 
some confidence in the use of isotropic TOV equations 
which assume P± = Pii . For the larger value of the mag- 
netic field studied, we find that the longitudinal pressure 
is still positive; however, the level of pressure anisotropy 
is quite large with P\\/P± — 0.4. In this second case, it is 
questionable whether isotropic TOV is sufficient to draw 
firm numerical conclusions; however, as we have demon- 
strated, the qualitative behavior and pattern of evolution 
observed with this larger magnitude of the magnetic field 
are similar to those obtained with the smaller magnitude 
magnetic field. 

Additionally, our results show clearly that the MIT bag 
model for a B 1 ^ 4 = 154 MeV obtained from an investi- 
gation of the adequate stability window [58| cannot re- 
produce the very massive neutron star recently detected 
[85|, not even if very intense magnetic fields are consid- 
ered. This is perhaps not surprising since most models 
which arc able to describe large mass stars rely on the ex- 
istence of additional repulsive (vector) interactions. That 
being said, it is interesting to consider the effect of a mag- 
netic field in the somewhat simple bag model in order to 
establish a baseline for other models. 

Finally, we note that although the results presented im- 
ply that quark stars with B c = 2 x 10 18 G are only mildly 
anisotropic, larger central magnetic fields cause the pres- 
sures to become increasingly anisotropic. Therefore, at 
such values of the magnetic field one should solve Ein- 
stein's equations in an axisymmetric metric which is de- 
termined self-consistently from the axisymmetric energy- 
momentum tensor for the star. Numerical solution for 
the axisymmetric case is needed and we are currently 
working towards this goal. Other issues to investigate 
in the future are, for example, the inclusion of the effect 
of the anomalous magnetic moment and to further refine 
the stability windows for the bag model in the presence 
of a background magnetic field. 
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